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Abstract 



We prove that the hnear term and quadratic nonhnear term entering a 
nonhnear elhptic equation of divergence type can be uniquely identified by 
' the Dirichlet to Neuman map. The unique identifiabihty is proved using the 

II , complex geometrical optics solutions and singular solutions. 
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■ 1 Introduction 



Let Q he a bounded domain in R" (n > 2) with a boundary dQ. We consider 
the fohowing nonlinear boundary value problem: 



(1.1) 



V ■C{x,Vu) = inn, 
u\dn = /• 



The nonlinear function C takes the following form: 

(1.2) C{x,q)=jq + Q{x,q):=jq + P{x,q) + R{x,q), x G fl, g G R" (or C"). 
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Here 7 G C^(0), and there exists Cq > such that 

(1.3) 7(x) > Co X G 
q) represents the quadratic nonhnearity: 

(1.4) P{x,q)=l cli{x)qkqi,--- , J2 

\l<k<l<n l<k<l<n 



Cki(x)qkqi 



We suppose that c^i G C^(0) are real valued functions and 

(1.5) hllcMH)' VkiWc^m <Ci {l<i<n,l<k<l<n) 

for some Ci > 0. R{x, q) G C^ljl x /f) with H := {q e R" or C" ; \q\ < h] {h > 0) 
represents the higher order nonhnearity, namely, it satisfies 



{ \d:,.R{x,q)\ <C2W^ 

(1.6) 



\d,^d,,R{x,q)\<C2\q\\ 



■3^ Ik' 

d,^dq,R{x,q)\<C2\f 
(1 < j,k < n,x £ Q,q e H) 

for some C2 > independent of x and 

The following theorem is probably well-known. However, since we are not able to 
find a proper reference, we give a short proof of it based on the contraction mapping 
principle in the Appendix at the end of this paper. For the elasticity equation, a 
proof using the implicit function theorem can be found in Q . 

Theorem 1.1 Let n < p < 00. There exist e and 5 < h/2 such that for any 
f G W^^^^^'P{dil.) satisfying ||/||^2~i/p,p(gf^) < e, j| j admits a unique solution u 
such that \\u\\iY2,p(^Q-^ < 6. Moreover there exists C > independent of f such that 

(1-7) lhllvi/2,P(C) < C\\f\\^y2-l/p,p(^QQy 

We define the Dirichlet-to-Neumann (DN) map Ac{f) for / with ||/||vi/2-i/p,p(gf^) 
< e to be 

(1.8) Ac(/) := C{x,Vu)\9n ■ v G W^-^'^^^m) 

where u is the unique solution of ( |1.1D such that ||tt||iy2,p{n) < ^■ 

In this paper we consider the inverse boundary value problem to identify C by 

means of the DN map Ac*. We are particularly interested in finding the linear term 

7 and the second order nonhnearity P{x,q). 

This inverse problem has interest in it's own right as the one to find conductivity 

distribution when the conductivity varies depending on the currents. It may also 
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considered as a simplified model for the nonlinear elasticity equation. In elasticity 
7(x) corresponds to the Lame moduli and c^;(2;) can be thought as the third order 
elasticity tensor. In acousto-elasticity, this higher elasticity tensor is important [^. 
We obtain the following unique identifiability theorem for 7 and P. 

Theorem 1.2 Suppose that n > 3. Let C^^\x,q) = 'jjQ + Q^^\ j = 1,2, where 
qU) = p(j){x,q)+R(^\x,q), 




P^^\x,q)={ c^}\x)qkqu--- , ^ 

l<fc<Kn 



and R^^^ satisfy ^1-^) , (1-5), and (L6), respectively, for the same constants Cq,Ci,C2 
and h. If 

(1.9) A^d) (/) = Ac.(2) (/) for all complex-valued small f G W'^-^/P'P{dn), 
then 

(1.10) 71(2;) = 72(2^) o.nd c^^i\x) = c^^i\x), a; G 
for 1 < j < n and 1 < k < I < n. 

The dimensional restriction n > 3 is imposed since we are using the complex 



geometrical optics solutions of Sylvester-Uhlmann |13| to prove the Theorem. 

In most models of nonlinear elasticity, the higher order tensors do not depend 
on X. So we consider the same inverse problem when the coefficients c^^ are con- 
stants. In this case we obtain the following uniqueness theorem including the two 
dimensions. One thing to be noted is the condition (|1.9| ) (and ( |1.11| ) below). Since 
the equation considered in this paper is nonlinear, A(^(i) (/) = A^(2) (/) for all real- 
valued small / G W'^~^^P'P{d^l) does not imply the same for complex-valued small 
/ G W'^~^^^'^{dO,). Therefore, in the situation where the data Ac(/) for only real- 



valued small / are available. Theorem 1.2 may not be applied. The following theorem 
uses only real-valued Dirichlet data. 



Theorem 1.3 Suppose that n>2. Let C^^\x,q) as in Theorem l.i except that the 



coefficients c^^^^ are constants. If 

(1.11) A^(i) (/) = A^(2) (/) for all real-valued small f G W'^-^/P'P{dn), 
then 

(1.12) ^i(x)=72(rE) and c«^' = 4f 
for 1 < j < n and 1 < k < I < n. 
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There have been some related works on the identification of nonlinear terms en- 
tering partial differential equations. Sun proved the global uniqueness for the non- 
linear conductivity 7(3;, u) when the Dirichlet data are complex valued and small 
|11|. Nakamura and Sun considered the similar problem for nonlinear elasticity 
model of St.Venant-Kirchhoff jl^. They proved that even in the presence of non- 
linearity, the linear term can be identified uniquely by means of DN map. Sun and 
Uhlmann proved the global uniqueness up to diffeomorphism fixing the boundary for 
the two dimensional nonlinear anisotropic conductivity A{x,u), and also the same 
result for the 3 dimensional analytic conductivity [|l2| . Quite recently, Isakov proved 
the unique identifiability of the nonlinearity c{u,p) entering the quasilinear elliptic 
equation Au + c(n, Vn) = (and corresponding parabolic equation) by means of 
DN map on dCl Q. G. Uhlmann informed us that Hervas and Sun proved a result 
related to ours in the two dimensional case ||6|. They proved that the constant co- 
efficients nonlinear terms with extra symmetry can be identified from the DN map 
defined for complex valued small Dirichlet data. 

Theorem and Theorem 1.2 are proved by investigating the first and second 
terms in the asymptotic expansion of Ac near 0. The first term is A^, the DN map 
corresponding to the conductivity 7. This fact was also observed in [10| and 0. 
We derive certain cubic relation from the second term. This is included in Section 
2. Then using the complex geometrical optics solutions of Sylvester-Uhlmann |13|, 
we prove Theorem 1.2 in Section 3. We prove Theorem 1.3 using the singular 
solutions of Alessandrini |2|. Since the coefficients to be determined are constants, 
boundary determination is sufficient. This is the reason why the singular solutions 
are effectively used. The proof is included in Section 4. Appendix is to prove 



Theorem 1.1 



2 Asymptotics of DN map 

Throughout this paper || W^^p denotes the VF'^'P-norm on Q, and || ||p = || ||o,p. Also, 
C > denotes the general constant in estimate independent of the functions being 
estimated. We will use the following estimates repeatedly: if k > 1, p > n, and 
n = 2, 3, then 



(2-1) II^^^^IU.P < C'lklU.plbllfc, 



This inequality holds by the Sobolev embedding theorem (see [^). We also use the 
chain rule for the derivative of the composition R[x, Vu{x)) for u E W'^'^{Q,) and its 
estimate given in Chapter II, section 3 of |14]. 

Suppose that / G W'^~^/P'P{dQ), p > n, and let t G be a small parameter. 
Let u*^*^ be the solution of 

(V ■C{x,Vu) = inQ, 
1 u\an = tf 
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such that ||u^*^||2,p < 6. Then there exists a constant C > such that 

(2.2) |kW||2,p<Ct||/|| 

Let ui = ui{f) £ W'^'^i^l) be the solution of 

fv • (tV-ui) =0 in 17, 
\ui\dn = /• 

Then by the regularity of the Dirichlet problem for elliptic equations, there exist a 
constant C > such that we have 

(2-3) ||mi(/)||2,p < C'll/llvF2-vp,P(an) 

for any / G W'^~^^P'P{d^l). Let U2 be the solution of 



V • (7Vn2) = -V • P{x, Vui) in n, 
U2\dn = 0. 



(2.4) 



It then follows from the regularity of elliptic equations, ( |2.1| ), and (^) that 

\\u2hp < C||V-P(x,Vni)||p 
<C\\Vu,\\l^ 

We now define v*-*^ by 

(2.5) -uW =t(ui + tt;W). 

Then, we have from and (HJ) 



t^v(^^\\2,p < ||nW||2,p + t||ni||2,p 
Lemma 2.1 There exist to and Cf depending on f, not on t, such that for all t < to, 

(2.6) Wv'^^^ -U2\\2,p<Cft. 

In particular, we have 

(2.7) ||^^(*)||2,p<C/. 
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Proof. Put w^^^ := v*-*^ — U2. Then a straight-forward computation shows that 

Jv • 7Vt(;(*) = -tgW in 0, 

where 

= t-iy • \p{x, Vni + tVw(*)) - P{x, Vni)j + t^^V • R{x, tVui + t^Vt;^). 
It follows from (|L§, (|lD, and 1^ that 

115^*) lip < C(||ni |b,p|bW ||2,p + t||^^(*) lllp + ||ni lllp + t'\\v^'^ |||,). 
It then follows from (^) that 

||^(*)-«2||2,p = lk(*)||2,p <Ct\\g(')\\p 

<C{t\\v(')\\2,p + t). 

Thus, if to is so small that Cto < |, then by (^) we obtain 

lb^*)||2,p<2(||tX2||2,p + l)<a 

Then ( ^.61) follows from This completes the proof. □ 



Lemma 2.2 For a given f G W^-'^/P'P{dn), we have 

(1) lim-^[Ac{tf)-tA,{f)] = 0, 

(2) lim i[Ac(t/) - tA^if)] = V ■ [7Vn2 + P(x, Vtxi)]. 
The convergence is in the topology ofW'^-'^/P'P{dn). 

Proof. From (^]^) we have 

Acitf) - tA^if) = V ■ [Cix, V^^W) - tTVni] 

= u ■ [t^-fVv'^^^ + Q{x,Vu^^y)]. 
It follows from the trace theorem Q, ( ^.2| ) and ( |2.7D that 

||t-i(Ac(t/) -tA7(/))lli^i-i/p,p(af7) < Wt^Vv^'^ +t-'Q{x,Vu^'^)\\i,p 

< Ct. 

This proves (1). 
Since 

Q(a;, VnW) = t^P{x,Vui) + t'^0{\Vui\\Vv^''\ + |V7;W|2 + |Vni|3 + |Vt;W|3), 
(2) follows from This completes the proof. □ 
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3 Proof of Theorem 



Lemma 2.2 says that we can recover and u ■ ['j\/u2 + P{x, Vui)] on 90 from Ac. 

If A^(i) = Aq(2), then A^^ = A^^. It then follows from well-known results (|^, 
g, I, @) that 71 = 72 in Q. Put 7 = 71 = 72. By Lemma U (2), we have 

(3.1) iy-[-/Vu'^'^ +P'^^\xyui)]=iy[-fVu^^^ + P^-^\x,Vui)] on dQ 

where ui G W'^'^{Q) is a solution of V • (7Viii) = in 17, and is the solution of 
(13) when P = P(J). 

Let V G Ty^'^(O) be a solution of V • {"yVv) = in 17 with v\qq = g. Since 

(7) 

'^2 \d^i = 0, it follows from the divergence theorem that 

V ■ [7V4^'^ + P^^\x, Vui)]gda = [ [7V4^'^ + P^^\x, Vui)] ■ Vvdx 

Jn 

P^^\x,Vui) ■ Vvdx. 



an 



We thus have 



(3.2) / P^^\x,Vui) - Vvdx = / P^^\x,Vui) - Vvdx, 

Jn Jn 

in other words, 



f3 3) V V A'^'(x)^^—dx-Y^ V A^>Ux) — ——dx 

1=1 l<j<l<n ■' i=l l<j<l<n ■' 



for all ui,v G VF^'^'(O), solutions of V • 7VU = in 17. Therefore, it suffices to prove 
that if 

1=1 l<3<l<n ■' 

for all ui, w G iy^'P(17), solutions of V • ^Vu = 0, then 

(3.5) cii = Q, i = l,--- ,n, l<j <l<n. 

Let W2 G Ty^'^(17) be another solution of V • jVu = 0. Then ui + ti2 is also a 
solution of the equation, and hence we obtain by polarizing ( |3.4D that 

(3.0) E Ej„c.,(.,(|i||.|||^)£..^o 

1=1 l<J<(<n ^ J J ^ 
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for all ui,U2,v G VF^'P(O), solutions of V • jVu = 0. The arguments we made so far 
are true for n > 2. 

Let n = 3 and use complex geometrical optics solutions. Let k E R'^ and choose 
unit vectors ^ and rj in so that 

(3.7) k-^ = k-rj = (,-T] = 0. 
Then choose t, s > so that 

(3.8) t^ = K + s\ 

4 

Define p^^\ p^"^^ G by 

(3.9) P^^'> =tr] + i{^+sO, P^^^ = -tr] + - sO 



Then by the fundamental work of Sylvester-Uhlmann [13| there exist solutions 
j = 1, 2, of V • 7Vtt = of the form 

(3.10) Uj{x) =-f-'^/^eP^'^-''{l + iPj{x,p^^^)), xGQ 
where ijjj satisfies 

(3.11) ll^illL°°{Q) < and W'^Tpjh-oin) < C 

for some C independent of t. We apply these solutions Uj to (^) and obtain 

^ / 7(x)-^c},(x)^(x)e^'=-x 

+{^y^da-^/\i + ^i) + pf\i + ^i) + di^i)x 

(^1/2^.^-1/2^1 + ^^2 + /of (1 + V'2) + 5,-V'2)' 

= 0. 

Here dj = d/dxj. Set Q := ij + and note that 

lim = lim = — 

s— >oo s s^oo s 



'J' 
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Therefore, by dividing the both sides of ( 3.12 ) by t, taking the hmit s —>■ oo, and 
using (|3.11|) 



I, we get 



(3.13) /EE C,Cic)i{xMx)-'^{x)e^'-^dx = 0. 

•'^ i=l l<j<l<3 * 



1<J 

bmce holds for ah A; G R^, we have 

3 



(3.14) 



E 

i=l 



E oc/c^^ 

l<j<'<3 



dv 

-— (x) = 0, X eft. 

OXi 



We need the fohowing lemma the proof of which will be given at the end of this 
section. 

Lemma 3.1 Suppose n > 2. There exist solutions Vj G W'^'^{Q), j = 1, ■ ■ ■ ,n, of 
V • j'Vvj = in ft such that 



(3.15) 

Lemma ^]l] and ( |3.14| ) yield 



-g-^{x) ) 7^ 0, for almost all x £ Q. 



(3.16) 



E CjC/c;7(^)=0' xGn, i = 1,2,3. 



Note that ( gl6D holds for ah C in the set V := {C e C^j C • C = 0, |C| = V2}. If 
C = (0, z, V^z) for z G C with \z\ = I, then (eV. With this C, (HH) becomes 



(3.17) 



42 (^) - 43 (^) + V^cis ix) = o, x£n. 



Since c*; are real, we have 022(2^) = c|3(x) and C23 = 0. By substituting C 
{z, 0, \/—lz) and ^ = (z, \/—lz, 0) in order into ( |3.16| ), we obtain that 

(3.18) 4=42 = 43 and 4, = (j / /). 

Because of ( 3.18 ), ( |3.6D now takes the form 



(3.19) 



3 3 

EE 

i=i j=i 



avians 



Q ^ ^ f?Xj dxj dxi 



for all uu2,v. By the same argument as above (using U2 and i'), we can conclude 
that 



(3.20) 



E 



dxj 



ix) = 
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for all ui and C ^ 1^- It then follows from Lemma 3.1 that 

3 



(3.21) 



CjY.Qcij{x) = 0, j = 1,2,3, 



i=l 



and we conclude that Cjj = 0, i,j = 1,2,3. The proof is complete. 



□ 



Proof of Lemma p.lj We assume that n = 3. The same proof works for the two 
dimensional case since complex geometrical optics solutions exist in two dimensions. 
Choose C^^^ G V so that they are linearly independent over C, e.g., ^^^^ = 

(1,V^,0), C(2) = (1^0,^/^), C^') = (0,l,\/^). Let 

^,-(x) = 7-'/V^*''-"(i + V',(x,tC^^'))), xen, j = l,2,3 

as before. Then, we have 

'vvi\ /c«(i+vi)+o(t-^r 



It then follows from ( |3.11 ) that if t is sufficiently large, then 

det Vv2{x) I ^ 0. 

\VV3{X) 



This finishes the proof. 



□ 



4 Proof of Theorem [TT3| 

We need to prove that if 

1=1 l<k<l<n 

for all u,v £ M^^'^(O) real solutions of V • 'yVu = 0, then 

(4.2) cii = 0, i = l,--- ,n, l<k<l<n. 

Since holds only for real-valued solutions u, v, we can not use the complex 
geometrical optics solutions. However, since the coefficients c^^ are assumed to be 
constants, we can use instead the singular solutions of Alessandrini. We only deal 
with the three dimensional case. Two dimensional case can be proved in the same 
way. 
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Fix j and s < t. We will show that c^f = 0. Let e^, k = 1,2, 3, is the standard 
basis of R'^. Let a > 0, /3 > 0, + = 1. Define an orthonormal frame N,Ti,T2 
hy N = acs + (3et, T2 = (3es — aet, and Ti = where k ^ s,t s ^ t. In the case 
s = t, let = Cs, and Ti and T2 be and e/ where k ^ I ^ s ^ k. For x G 50, 
denote the unit outward normal to dfl at x. Choose xq S 50 so that ;/(a;o) = N. 
By translation if necessary, we may assume that xq = 0. 

Choose a solution f G C^(0) of V • 7V?; = so that div{0) 
can be proved easily that such a function exists. Then 



5ij (i = 1,2,3). It 



(4.3) divix) = 6ij + Oi\x\), x^O. 

Let U be an open neighborhood of 0. Then we have from (|4.1|) and ([1.3|) that 



(4.4) 



4; / dkudiudx 
i<k<i<n •'^^u 



< C I \x\\Vu\^dx 



for some C > 0. 

By PI, for every e > 0, there exists a solution u £ C^(0) such that 



u{x) 



1 



\x - eN\ 



+ zu(x) 



where 
(4.5) 



|Vu'(x)| < C\x-eN\ 



-2+a 



for some q > (a < 1) and C independent of e. In two dimensions we can use 
logjx — eA^|. Put re(x) := \x — eA^|~^ to make notations short. Substituting this 
solution u to ([4.4D, we get 



(4.6) 



V4/ / dkTATedx <C I \x-€N\ 
f,^i Jnnu Jo. 



-4+a 



dx. 



For convenience, put y := x — eN . From (|4.6|), we have 



(4.7) 



Y.<il fi^- [ \- - ^N^'-'dx. 
Jnnu \y\ Jn 



Suppose dQ near is given by 

aon[/ = { 6^1 + 6^2-9^(6,6)^^ I 161' + 161' < '5' } 
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where tp G C^, 95(0) = 0, and Vv?(0) = 0, for a fixed small 5 > 0. Let D : = 
{ ^iTi + 6^2 - ^sN I iCiP + |6P < 5^ < 6 < (J } and define <!> : D ^QnU hy 

^^iTi + ^2^2 - ^sN) := eiTi + 6T2 - (6 + '/'(a, 

By shrinking U if necessary, we may assume that ^{D) = Q CiU. Then we have 



after a coordinate change, (4.7) takes the form 



(4.8) 



k<l 



where z = ^{S,) — eN, ^ = S^iTi + ^2^2 — Cs-^ £ -C^- We now scale r/ — > ^ = e?7 to 
obtain 



(4.9) 



Jd J\v 



5/e 



CkCl 
\C\' 



dr] 



where ( = rjiTi + 772^2 — (??3 + e ^(p{er]i, er/2) + 1)A^. It is easy to prove, by scaling, 
that 

(4.10) lime [ |x - eiVr^+"dx = 0. 

Jn 

Let rj' := {r]i,r]2). Since 



one can see that 



e~^\ip{er]')\ < Ce\7]'\^ < C6\r]'\, for |r?'| < 6/e, 



CkCl ^ ^ 1 



|C|6 - |r?|4 + l 

for some C independent of e. Therefore, we obtain from the dominated convergence 
theorem 

(4-11) lim/ r%^d,= [ 

where w; = 771 Ti + 772 T2 — (773 + 1)A^ and Ri]_ is the upper half space. Let 

It follows from (|^, (gl^) and ( plj) that 

(4.12) E« = 0. 
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Suppose that s = t. In this case was chosen to be = e^. Thus Wk = rjk if 
k ^ s and Ws = — (r/^ + 1). Therefore easy computations show that 

if A: / /, 



^ if k = I s, 

— if k = I = s. 
I 2 



Thus by (p^ ) we get 
(4.13) 



24s + ^4k = 0, 5 = 1,2,3. 



Hence we conclude that 

(4.14) 4^ = 0, j,fe = l,2,3, 
and ( [4.12 ) takes the form 

(4.15) E« = 0- 

k<l 

Suppose that s < t. If k ^ s,t, then Ti = and hence Wk = Ws = f3rj2—a{r]3+l), 
and = —ar]2 — /3(r/3 + 1). Again by simple computations one can see that 



jst 
^kl 







a/3 



TT 



if (A:,/) / (s,t), k<l, 
if (A;,/) = (s,t). 



It then follows from ( [4.15 ) that for s <t 

(4.16) ci^ = Q, j = 1.2,2,. 

By ( [4.14 ) and ( 4.16] ), the proof is complete. 



□ 



5 Appendix- Proof of Theorem 

Replace R{x,q) by R{x,q) := x{q)R{x,q) + (1 - x{q))\<l? with x S C^(R" or C'') 
satisfying < xil) ^ 1) x('?) = 1 (kl ^ V^)) i\q\ > h). Here, for the complex 
vector q S C"" is identified with (Reg, luiq) G R^" and xil) is considered as a 
function of variables Reg, Img. Then, (|1.6| ) holds for any g S R" or C" with some 
new C2 determined by h, xil) and the old C2. Suppose that ||/||i4/2-i/p,p(sn) < e (e 
to be determined later). Let uq be the solution of 

^^^^ fV-7Vuo = 0, inJ7, 
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V --fVv + V ■Q{x,Vuo + Vv) = 0, in 0, 
v\dn = 0. 



If n is a solution of ( |1.1D , then v defined hy u = uq + v satisfies 
(5.2) 

Let be the solution operator for the problem 
(5.3) 



V ■^Vu = g e LP{n), in 17, 
u\dQ = 0, 



namely, L^^g is the solution of (^.3|). Then L^^ is linear and there exists C > 
independent of g such that 

\\L::^'^g\\2,p < c\\g\\p. 

Let X := {if ^ W'^'P{Q,) \ (p\dn = 0, ||93||2,p < 5 } (5 to be determined later) and 
define an operator A on X by 

A{ip){x) := L-^V ■ Q{x, VuQ + Vip). 

Then, it follows from ( |1.4| ) and ( |1.6D that 

\\A{ip)\\2,p <C||Q(Vno + V<^)||i,p 

(5.4) < C(||no||L + IMlp + Ikolllp + Mllp) 

<Cie^ + 5^ + €^ + 6^), 

and 
(5.5) 

\\A{ipi) - A{ip2)\\2,p < C\\Q{x,Vuo + Vipi) - Q{x,Vuo + Vip2)\\i,p 

< C{\\uo + Oifi + (1 - e)ip2\\2,p + II^^O + 0^1 + (1 - 0)ip2\\lp) 
\\V{(pi - (fi2)\\l,p 

<C{e + S + e^ + 6^)yi-ip2\\2,p 

for some < 6 < I. Therefore, if C{e^ + 6^ + + 6^) < 6 and C{e + 5 + e^ + 6^) < 1, 
then A is a contraction on X. 

By the contraction mapping principle, there exists a unique fixed point of A, say 
V. Then v is the solution of ( |5.3D . The estimate ( |1.7D follows from the equation and 
regularity of the linear elliptic equation. 

Finally, we take e and 6 more smaller if necessary so that they satisfy C^e + 5 < 
h/2. Then, u = uq + v is the solution of (1.1). 
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